ABSTRACT
A novel severe acute respiratory syndrome (SARS)-like coronavirus (SARS-CoV-2) has emerged as a human pathogen, causing global pandemic and resulting in over 400 000 deaths worldwide. The surface spike protein of SARS-CoV-2 mediates the process of coronavirus entry into human cells by binding angiotensin-converting enzyme 2 (ACE2). Due to the critical role in viral-host interaction and the exposure of spike protein, it has been a focus of most vaccines' developments. However, the structural and biochemical studies of the spike protein are challenging because it is thermodynamically metastable. Here, we develop a new pipeline that automatically identifies mutants that thermodynamically stabilize the spike protein. Our pipeline integrates bioinformatics analysis of conserved residues, motion dynamics from molecular dynamics simulations, and other structural analysis to identify residues that significantly contribute to the thermodynamic stability of the spike protein. We then utilize our previously developed protein design tool, Eris, to predict thermodynamically stabilizing mutations in proteins. We validate the ability of our pipeline to identify protein stabilization mutants through known prefusion spike protein mutants. We finally utilize the pipeline to identify new prefusion spike protein stabilization mutants.